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ABSTRACT 

Models for the formation of galaxies and clusters of galaxies require strong feedback 
in order to explain the observed properties of these systems. We investigate whether 
such feedback has observational consequences for the intergalactic medium, as probed 
in absorption towards background quasars. A typical quasar sight-line intersects one 
proto-cluster per unit redshift, and significant feedback from forming galaxies or AGN, 
heating the proto-cluster gas, will result in a large clearing of reduced absorption 
in the Lyct-forest. Such a gap could be detected at redshift <; 3 when the mean 
opacity is high. Feedback from Lyman-break galaxies in proto-clusters can be probed 
by the absorption lines produced in their winds. Strong feedback from galaxies has a 
major impact on the number and properties of absorption lines with column densities 
Am ~ 10 16 cm~ 2 . This feedback can be probed with multiple sight-lines and by 
studying the unsaturated higher-order lines of the Lyman series. Galactic winds from 
dwarf galaxies should break-up into clouds, in order not to over produce the number 
of absorption lines. These clouds can then coast to large distances. 

Key words: cosmology: theory - intergalactic medium - hydrodynamics - large-scale 
structure of universe - quasars: absorption lines - galaxies: formation 



1 INTRODUCTION 

In cold dark matter (CDM) dominated cosmogonies struc- 
tures grow by gravitational amplification of small primor- 
dial density perturbations. High-density regions decouple 
from the Hubble expansion and form gravitationally bound 
systems with a range of masses up to a characteristic cut- 
off mass that increases with decreasing redshift (Press & 
Schechter 1974). Early on, pressure forces are unable to sep- 
arate the gas from the dark matter. Galaxies form when 
baryons that are dragged into the dark matter potential 
wells are shocked to sufficiently high temperatures and den- 
sities that the gas can cool radiatively (Rees & Ostriker 1977, 
White & Rees 1978). 

This hierarchical picture of galaxy formation predicts 
the formation of too many small systems and therefore con- 
flicts with observations of the faint-end slope of the galaxy 
luminosity function. In order to resolve this discrepancy, 
feedback is assumed to suppress the formation of most small 
galaxies (Larson 1974, Rees & Ostriker 1977, White & Rees 
1978). The required feedback could be energy input by su- 
pernovae (SNe) associated with star formation (e.g. Dekel & 
Silk 1986). In more massive galaxies, feedback is assumed to 
be able to prevent the catastrophic loss of angular momen- 
tum of forming discs, which causes discs in simulations to be 



much smaller than observed discs (e.g. Navarro & Steinmetz 
1997, Weil, Eke & Efstathiou 1998). 

On even larger scales, the hierarchical picture based 
solely on gas shocking and radiative cooling also seems to 
conflict with observations of groups and clusters of galaxies. 
Without some type of feedback, shocked gas in the halos of 
galaxy groups will over produce the X-ray background (Wu, 
Fabian & Nulsen 1999b, Pen 1999). Feedback is also invoked 
to explain the observed deviation from self-similarity of the 
X-ray properties of hot gas in galaxy clusters (Kaiser 1986, 
1991). 

Therefore, both theoretical and observational argu- 
ments suggest that forming halos are not just passively ac- 
creting material from the intergalactic medium (IGM), but 
may be a significant source of energy input into their sur- 
roundings. Is this energy injection a fairly local process, only 
affecting gas within a small fraction of the virial radius of 
the halo say, or do galaxies stir-up a significant fraction of 
the IGM? 

This question can be approached from a different an- 
gle, by studying the properties of the IGM in absorption to- 
wards background quasars. Neutral hydrogen along the line 
of sight produces the Lya absorption lines observed as the 
quasar's 'Lya' forest (Bahcall & Salpeter 1965; Gunn & Pe- 
terson 1965, see Rauch 1998 for a review). Several theoretical 
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models examined whether collapsed objects could produce 
the observed absorption until it was realized that the IGM 
itself gave rise to the majority of the lines (McGill 1990, Bi, 
Boerner & Chu 1991). Numerical simulations of CDM dom- 
inated cosmologies in which the IGM is photo-ionized by 
the UV-background from quasars, have been very successful 
in reproducing many observed properties of the Lya-forest 
(Cen et al. 1994, Zhang, Anninos & Norman 1995, Hern- 
quist et al 1996, Theuns et al 1998b; see e.g. Efstathiou, 
Schaye & Theuns 2000 for a review). For example, Theuns 
et al (1998a) showed that this model naturally reproduces 
the evolution in the line densities per unit redshift and unit 
column density, from redshift to 4. The superb quality 
of the data that this comparison is based on increases our 
confidence that the current paradigm of the Lya-forest - 
absorption by neutral hydrogen in the intervening IGM - is 
basically correct. 

In these simulations of the Lya-forest that do not in- 
clude any feedback from forming halos (other than the pres- 
ence of a UV-background), strong Lya absorption lines (col- 
umn densities > 10 15 cm~ 2 ) are invariably produced close 
to forming halos. If feedback were indeed important, as ex- 
pected from the earlier theoretical arguments and from the 
actual observations of strong galactic winds associated with 
star bursts (e.g. Heckman, Armus & Miley 1990), one would 
expect the properties of these stronger lines to be at least 
partly influenced by the feedback process. 

If galaxy feedback is required to explain the observed X- 
ray properties of groups and clusters, then real proto-clusters 
might be very different from the proto-clusters in the simula- 
tions that do not include feedback. Consequently, one might 
expect the absorption properties of the IGM where the line 
of sight intersects a proto-cluster to be rather different in 
the presence of feedback. 

Metal lines (e.g. Civ) associated with Am ~ 
10 14 ' 5 cm" 2 column density lines have been detected in high 
signal-to-noise QSO spectra (e.g. Cowie et al. 1995). Even if 
the lower-density gas were polluted by metals at the same 
level, the associated Civ metal lines would generally be too 
weak to be detected. Ellison et al. (2000) searched for Civ 
in a very high quality spectrum on a pixel by pixel basis 
and were able to show that there must be many more metal 
lines than can be detected directly. Schaye et al. (2000a) 
performed a similar search for Ovi and found evidence for 
oxygen enrichment down to TLyc* ^ 10 _1 atz~ 2-3, which 
corresponds to underdense gas. This is again circumstantial 
evidence that a galaxy can influence the gas around it out to 
a significant distance. The case is not so clear cut, because 
the metal pollution could be the result of an earlier epoch of 
population III star-formation. But if galactic winds were the 
source of the pollution, then these winds might change the 
properties of the low density IGM as observed in absorption. 

The aim of this paper is to investigate to what extent 
the current picture of the Lya-forest can incorporate a sig- 
nificant amount of feedback, without destroying the good 
agreement with observations. In what follows we will also 
use numerical simulations of the Lya-forest. We simulate a 
vacuum energy dominated, flat cold dark matter model, with 
density parameters (fi m ,S7A) = (0.3,0.7) (e.g. Efstathiou 
1999 and references therein), Hubble constant h — 0.65 
(Freedman et al. 1999) and baryon density Q, b h 2 = 0.019 
(Buries & Tytler 1998) of which a fraction Y = 0.24 by mass 



is helium, normalized to the abundance of clusters (Eke, 
Cole & Frenk 1996). We assume the IGM to be photo-ionized 
and photo-heated and chose the spectrum of the ionizing 
sources such that the simulation reproduces the thermal evo- 
lution and mean absorption as determined from observations 
by Schaye et al. (2000b; the gas temperature as function of 
redshift follows the dashed line in Figure 6 of that paper). 
The cubic simulation box is 20 Mpc on a side (co-moving) 
and contains 256 3 particles of each species (i.e. dark matter 
and gas particles; note that the designer model in Schaye 
ct al. (2000b) was a smaller simulation but with the same 
thermal history), providing sufficient resolution to simulate 
the gas reliably (Theuns et al. 1998b). Without feedback, 
this simulation reproduces most of the observed properties 
of the Lya-forest, for this choice of currently popular values 
of the cosmological parameters, it includes star-formation, 
whereby gas at a density contrast p > 80 (p) is converted to 
'stars' , when the temperature drops below 2 x 10 4 K. No 
feedback (other then a UV-background) is included. 

This paper consists of two parts. In the first part we 
assume that the gas in a proto-cluster is pre-heated at some 
early time, and discuss the corresponding observational ab- 
sorption signature if the line of sight to a quasar intersects 
the proto-cluster. In the second part, we consider feedback 
effects from smaller systems and examine their impact on 
the IGM. 



2 FEEDBACK ON CLUSTER SCALES 

We begin this section by briefly reviewing the arguments 
that suggest that the entropy of proto-cluster gas was in- 
creased significantly by some non-gravitational process at 
an early time Zi n j > 2. We then argue that this gas is too 
hot to be confined by a galactic potential well and hence will 
fill a considerably fraction of the volume of the proto-cluster. 
The probability that a quasar sight-line intersects such a hot 
bubble is ~ 1 per unit redshift, and it will appear as a re- 
gion of decreased absorption in the Lya-forest. This should 
be relatively prominent especially at high redshift z > 4 
when the mean Lya-opacity is high. Other processes, such 
as proximity effects from AGN and/or galaxies in the proto- 
cluster, could contribute to clearing a gap in the absorption. 
Observations of metal absorption lines can be used to exam- 
ine the thermal properties of the medium, and hence distin- 
guish a proto-cluster gap from a void, and make inferences 
about the nature of the sources producing the feedback. If 
the proto-cluster gas is in a two-phase medium, with cold 
clouds pressure-confined in the hot phase, then the gaps may 
not be completely empty. 



2.1 General considerations 

Clusters of galaxies are the most massive virialised systems 
in the universe, with typical total masses of about 1O 15 M0, 
of which 10-30 per cent is baryonic. The latter component 
is predominantly in the form of a hot X-ray emitting intra- 
cluster medium (ICM), with temperature and X-ray lumi- 
nosities typically of order T ~ 10 s K and Lx ~ 10 45 ergs -1 , 
respectively. 

The intracluster gas is believed to be shock-heated as it 
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falls into the cluster potential well. If the initial gas temper- 
ature is negligible, then a relation Lx oc T 2 is predicted 
from simple scaling arguments (Kaiser 1986), which dis- 
agrees with the observed result, Lx oc T a , with a = 2.6- 
3 (Jones & Forman 1984, David, Forman & Jones, 1991, 
Lloyd-Davies, Ponman & Cannon 2000). Consequently, it 
was suggested that the ICM gas is 'pre-heated' by some non- 
gravitational process (Kaiser 1991; Evrard & Henry 1991; 
Navarro, Frenk & White 1995; Cavaliere, Menci & Tozzi 
1997, Ponman, Cannon & Navarro 1999). The existence of 
such an entropy 'floor' breaks the self-similarity and helps 
to explain the Ly^-T relation (Kaiser 1991). Since an initial 
entropy is expected to have a larger effect in clusters with 
lower temperature, the level of pre-heating may be revealed 
by measuring the specific entropy of gas in X-ray clusters as 
a function of cluster temperature. The effect of pre-heating 
would be largest in groups of galaxies, which would emit 
far less X-rays when pre-heated since the hotter gas is less 
dense, thereby avoiding the problem of over producing the 
X-ray background (Wu, Fabian & Nulsen 1999b, Pen 1999). 
Defining a specific 'entropy' 



T 

7273 



(1) 



where T is the temperature of the gas and n c is the electron 
number density, then the level of pre-heating inferred from 
the observations is on the order of 



so 



lOO/i" 



-1/3 



keVcm 



(2) 



(Ponman, Cannon & Navarro 1999, Lloyd-Davies, Ponman 
& Cannon 2000, Balogh, Babul & Patton 1999, Wu, Fabian 
& Nulsen 2000). We will use sioo = s /(100/i" 1/3 keVcm 2 ). 

At present it is not yet clear which sources dominate the 
heating of the ICM. The energy required to reach a given 
entropy floor increases with the gas density as E oc n 2 / 3 , 
which makes it easier energetically to heat the gas before it 
collapses into a potential well. If supernova explosions and 
stellar winds, e.g. associated with episodes of rapid star- 
formation, are the main sources of energy, then pre-heating 
is likely to have occured relatively early on. Indeed, recent 
observations show that large numbers of stars are being 
formed at z ^ 3 (Steidel et al. 1996) in strongly clustered 
galaxies (e.g. Giavalisco et al. 1998, Adelberger et al. 1998) 
which are likely progenitors of the bright galaxies in present- 
day clusters (Mo & Fukugita 1996, Baugh et al. 1998, Wech- 
sler et al. 1998, Mo, Mao & White 1999). Semi-analytical 
models for the formation of cluster elliptical galaxies also 
predict that most of their stars formed at z <; 2 (Kauff- 
mann, White & Guiderdoni 1993, Kauffmann et. al. 1999). 
In addition, Mushotzky & Loewenstein (1997) argued that 
the observed lack of evolution in the iron abundance of clus- 
ter gas out to z ~ 1 also suggests that enrichment occured 
early. If, on the other hand, energy input from AGN is the 
main energy source, then pre-heating is also likely to have 
occured at z ^ 2, because the number density of QSOs peaks 
around z — 2-3 (Shaver et al. 1996). Therefore, for both 
types of potential energy sources, stars or AGN, we expect 
that the major part of pre-heating took place before z ~ 2. 

The inferred entropy floor Eq. (^|) implies that the ICM 
was heated to some high temperature at an early time. At 
a redshift z, the number density of electrons, n e , of fully 



ionized proto-cluster gas of primordial abundance (with a 
helium fraction of 24 per cent by mass) , at over density 6 is 



1.2 x 10 _5 (1 + 5) 



0.019 



In order to achieve the required entropy 
region must be heated to a temperature 



(3) 



the gas in the 



T 



0.06 



SlOO 



(1 + 5) 



2/3 



n b h 2 x 2/3 



0.019 



(^)Uf"-< 



(4) 



where f s and <j) are the heated mass and volume fractions 
of gas, respectively. The implied temperature is quite high 
(T £ 10 6 K; lkeV^ 1.2 x 10 7 K) if pre-heating occured at 
z > 2, even if the heating occured close to the mean density 
5 ~ 0. At z < 1, the required T is lower (T ~ 10 5 K) if 
5 = 0. However, for a cluster to form at the present time, 
the redshift at which it turned around is about 1, in which 
case the gas density when the pre-heating took place must 
have been 5 ~ 10 over the volume of the proto-cluster, so 
we again obtain T £ 10 6 K. 

If pre-heating occured inside galaxy halos, then /„ and 
cf> are small and the required temperature Eq. (Q) is higher 
than the virial temperature of the halos. Such hot gas cannot 
be confined to a galaxy halo and therefore, irrespective of the 
heat source, the pre-heated gas must have quite an extended 
distribution, filling a large fraction of the region occupied by 
a proto-cluster. 

If pre-heating is due to SNe, one can estimate the to- 
tal amount of stars needed to heat the proto-cluster to the 
required temperature. Suppose the gas mass of the proto- 
cluster is Mgas, and /*M gas is the mass of stars responsible 
for pre-heating. The total energy released by the SNe asso- 
ciated with those stars is 



E = 



-A£? 5 ilO ai erg, 



(5) 



where E51 is the energy output per SN in units of 10 51 erg 
and v is the stellar mass needed to produce one supernova 
(we will write v = lOOi^ioo Mq). If a fraction e of all of the SN 
energy is available to heat the gas, then the gas temperature 
will be 



kT w 2.1e 



^100 



keV. 



(6) 



It then follows that /* must be about O.le -1 to 0.2e~ in 
order to heat the proto-cluster gas to a few million degrees. 
This fraction is consistent with the observed fraction of bary- 
onic mass in stars (~ 10 per cent) and the metallicity of X- 
ray emitting gas in clusters, if the efficiency is high, e ~ 1. 
However, if only a small part of the feedback energy goes 
into heating the gas, then supernova explosions and stellar 
winds may not be sufficient to heat the gas to the required 
temperature, and heating from AGN may have to be invoked 
(e.g. Wu, Fabian & Nulsen 2000). 

Another requirement for pre-heating to work is that the 
heated gas does not lose the injected entropy radiatively. The 
time scale t coo \ = T/T for radiative cooling is 



1.0 x 10 1 



26A_ 23 



n b h 2 

0.019 
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yr, 



(7) 



where Tq is the gas temperature in 10 6 K and A(T) = A_23 x 



is the radiative cooling rate. For T> 



10 



~ 23 v / 3¥ergcm 3 s 



10 K, cooling is dominated by thermal Bremsstrahlung and 
A_23 ~ 1, whereas at lower temperatures > 10 4 K line- 
cooling becomes important and the rate can increase up 
to A_23(T) ~ 10 2 depending on the metallicity and the 
temperature of the gas. At high temperatures and for ten- 
uous gas, cooling will be dominated by Compton losses on 
the micro-wave background radiation, and the cooling time 
icooi ~ 1 x 10 10 [(l + z)/4]~ 4 yr. Comparing these times with 
the Hubble time 



tu 



3.4 x H^ft/O.^)- 1 (n m /0.3)~ 1/2 (j-t?-) yr , (8) 

(here and in what follows, we use the high-redshift approxi- 
mation H(z) fa Ho£lm 2 (l + z) 3 ^ 2 , valid for a cosmologically 
flat universe at z > 2) and using Eq. (Q), we see that the 
cooling time is equal to the age of the Universe when S ~ 100 
at z = 3. Again this shows that pre-heating should occur in 
gas at low over density, since otherwise the injected entropy 
will be lost radiatively. 



2.2 Observational signatures 

In the previous section we argued that pre-heating of gas in 
clusters is likely to have occured before z > 2, and that the 
gas temperature is fairly high, T Jj 10 6 K. We also argued 
that the hot gas should fill a relatively large fraction of the 
volume occupied by the proto-cluster, since the hot gas in 
unlikely to be confined by a galactic potential well. Because 
of the high temperature, the gas in a proto-cluster will be 
collisionally ionized thereby producing a region of reduced 
absorption in a QSO sight-line that intersects the cluster. 
In addition, proximity effects due to AGN and/or galaxies 
may further reduce the amount of absorption. What are the 
expected sizes and rates of incidence of such 'gaps' ? 

The rate of incidence of objects with co-moving density 
n and physical radius R is 



dN 

dz 



nnR 



H(z) 



(9) 



The physical radius of a proto-cluster of mass M = 
lO 15 Afisfc -1 M , which is at an over-density S at redshift 
2, is 



R « 10!T2~ 1/3 (1 + <5)~ 1/3 (1 + z^Ml^hT 1 Mpc, 



(10) 



and consequently the average number of times a sight-line 
will intersect such a a proto-cluster is (for z <; 2) 



^ « 1.9n 4 (fi m /0.3)- 7/6 (l + ( 5)- 2/3 M 2 5 /3 / ' + ' 



-3/2 



,(H) 



where n± is the co-moving number density of proto-clusters 
in units 4 x 10~ 6 h 3 Mpc~ 3 . For Abell clusters with mass 
A/15 0.5, n 4 ~ 1. 

The differential Hubble velocity V = H(z)R over a re- 
gion of this size is 



V 



1.6 x 10 3 (a„/o.3) 1/6 (i^ 1 ) 

(l + 5)- 1/3 M 1 1 5 /3 km a -1 , 



1/2 



and the corresponding wavelength range in a Lya spectrum 
is 



AA w 26 



(a „/o.3) 1 '' 6 (i±i) 3/2 (i + «5)- 1 



> 3 Mli 3 k. 



(13) 



(12) 



The observational signature of a proto-cluster will be smaller 
than this, since peculiar velocities will compress the cluster's 
extent in redshift space. The size of the cluster given by 
Eq. ([[3]) is large compared to the distance between lines at 
z ~ 3, and even more so at higher redshifts where the mean 
opacity in the forest is higher, as demonstrated in Fig. |l| 
This figure shows a simulated Lya-spectrum of a z cm = 
3 and a z cm = 4 quasar, extending over a velocity range 
typical for observed quasars. (Here and below, velocity v is 
related to wavelength A by dv/c = dA/A). Such spectra do 
contain 'gaps' with little absorption, which are due to large 
voids resulting from long- wavelength perturbations, but only 
very rarely are such voids of the size given by Eq. (|l2|), 
especially at the higher redshift where the mean absorption 
is higher. The simulation used to make Fig. |l| lacks large- 
scale modes due to it's finite size - 20 Mpc - and this may 
lead to an underestimate of the sizes of the largest voids. A 
more quantitative comparison against data could take such 
long-wavelength modes into account by using for example 
the log-normal model of Bi & Davidsen (1992). 

Thus, if pre-heating of the ICM occurred at z <; 3, its 
signature may be apparent in a QSO spectrum that inter- 
sects a proto-cluster. Therefore, QSO absorption line spectra 
are potentially powerful probes of feedback on cluster scales, 
since an average quasar sight-line will intersect a proto- 
cluster region of order ~ 1 times, according to Eq. (|ll|). 

However, there are complications in the model predic- 
tion. At low redshifts z < 1, proto-clusters may have already 
began to contract relative to the Hubble flow, and the clear- 
ings in the forest produced by the hot regions in QSO spec- 
tra are smaller than when assuming pure Hubble flow. Since 
the line density in the Lya-forest decreases with redshift, it 
becomes increasingly difficult to associate a clearing in the 
forest with a proto-cluster, since such gaps naturally appear 
as a consequence of the emptying of voids and the overall 
decrease in the Lya optical depth. 

Another complication is that some cold gas may exist 
in the hot proto-clusters and produce Lya absorption lines, 
so the clearings may not be empty. Cold gas^] is likely to 
co-exist with the hot phase in the form of pressure- confined 
clouds, in analogy with the multi-phase model of the in- 
terstellar medium (ISM, McKee & Ostriker 1977). In these 
models of the ISM, cold clouds are pressure-confined by a 
hot medium generated by supernova explosions. Such mod- 
els have also been considered for the Lya-forest itself (Os- 
triker & Ikeuchi 1983) as well as for the QSO absorption 
systems associated with galaxies (Mo 1994; Mo & Miralda- 
Escude 1996). 

We can estimate the properties of these cold clouds, 
by assuming pressure equilibrium between the hot and cold 
phases, 



* The sight-line may also intersect a galaxy. If the cluster contains 
10 3 galaxies, then the covering factor of intersecting one within 
30/i -1 kpc is 10 3 x [0.03(1 + z)/10] 2 ~ 0.1 for a R = lO/i^Mpe 
cluster at z = 3. 
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Figure 1. Simulated Lya spectra at redshift 2 = 3 (top) and z = 4 (bottom), with mean absorption scaled to that of observed spectra 
(Schaye et al 2000b). Vertical lines are drawn at intervals given by the expected size of a cluster, Eq. (|l2|). Intervals with little absorption 
over a region the size of a proto-cluster are very rare, especially at z = 4 when the mean opacity is high. A spectrum of this length 
contains of order one proto-cluster, which would stand-out as a region of decreased absorption if the proto-cluster gas were hot due to 
entropy injection. 



n c T c = n h T h . (14) 

If the cold clouds are heated by the general UV-background, 
their temperatures will be typically ~ 10 K. If f c = 
M c /Mg as is the fraction of gas in the cold phase, and m c 
is the typical mass of a cold cloud, then 



C=f e 




is the covering factor of the clouds. The mass of cold clouds 
in a hot medium is limited by various physical processes 
(evaporation, hydrodynamical and gravitational instabili- 
ties; see e.g. Mo & Miralda-Escude 1996), and typical values 
are m c ~ 1O 6 M0. Equation ( |i~5| ) shows the covering factor 
to be ~ 1 requires that about 10 per cent of the proto-cluster 
gas is in the cold phase. In that case, each time we observe 
a clearing in the forest, produced by a proto-cluster, one 
absorption line produced by the pressure-confined clouds is 
expected. If this component of absorption lines is indeed 
important, the proto-cluster gap in the spectrum would be 
smaller. 

The HI column density of a pressure-confined cloud is 
in the optically thin case 

where P is the pressure of the confining medium, R is the 
cloud's radius, T c is the temperature of the cold medium, 
and J21 is the ionization flux at the Lyman limit, in the 
usual units of 10 -21 ergs -1 cm -2 sr" 1 Hz -1 (see Mo & Mor- 
ris 1994). For typical values P ~ 20 (corresponding to 
T h = 10 6 K and n h ~ 2 x 10 -5 cm -3 ), the size of a cold 
cloud (with mass 1O 6 M0) is about 1 kpc, and the typical 
HI column density of such cloud is thus JVhi ~ 10 lj ' 5 cm -2 . 
Therefore these pressure-confined clouds can produce strong 
Lya absorption lines. Since small and diffuse clouds will 



evaporate in the hot environment, weak Lya lines are not 
expected in the hot gap. 

We will show below (Fig. ^) that the typical strong lines 
of column density Nm ~ 10 16 cm -2 produced in simulations 
of the Lya-forest are the result of the overlap in redshift 
space of several filaments. The Lya-profile of such a line is 
saturated but the substructure can still be probed by the 
higher-order lines (e.g. Ly/3). In addition, the extent of such 
absorbers perpendicular to the line of sight, which can be 
probed using double or lensed quasars, is very much larger 
than for the ~ lkpc pressure-confined clouds. Therefore it 
is possible to distinguish strong lines produced by a single 
cold cloud, as described above, from those that are typical 
in gravitational structure formation. 

If pre-heating is of stellar origin, then the ICM may 
be enriched in metals by supernova explosions and stellar 
winds. The mass- weighted mean metallicity of the proto- 
cluster gas is Z ~ yf t , where y is the metal yield per 
unit stellar mass. Since y ~ Zq, we have Z ~ f*Z®, 
i.e. the proto-cluster medium can be enriched to a level 
of 0.1-0.2 Zq, consistent with the observed metallicity of 
the present ICM (e.g. Renzini 1997). The absorbing clouds 
confined in the ICM are presumably enriched to a similar 
level. If this is the case, we would expect to observe strong 
metal lines associated with the Lya clouds. Using the typ- 
ical density and size of the cold clouds, and assuming the 
UV-background intensity has a power law form J(u) oc v~ x 
with amplitude given by J21 ~ 1, the associated CIV ab- 
sorption systems are expected to have column densities of 
10 13 -10 14 cm , comparable to the observed CIV column 
densities associated with the strong lines in the Lya-forest 
(e.g. Cowie et al. 1995). If the hot ICM is polluted by met- 
als as well, then a proto-cluster gap may be associated with 
very broad metal lines, due to the high gas temperature 
and differential Hubble expansion across the region. Such a 
characteristic signature may help in distinguishing a proto- 
cluster gap from a void. 

If, on the other hand, heating and metal enrichment 
have different origins (e.g., heating was done by AGN), met- 
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als produced by supernovae and stellar winds may be con- 
fined to small regions around individual galaxies. In this 
case, a large fraction of the heated ICM may remain at low 
metallicity, and the strong Lya lines produced by pressure- 
confined clouds may not be associated with any strong metal 
lines. Thus, by measuring the metallicities of the strong lines 
produced in the proto-cluster regions, we may hope to obtain 
constraints on the sources responsible for the pre-heating. 
Proto-cluster regions may be identified as large concentra- 
tions of galaxies in the galaxy distribution. Analyses of QSO 
absorption spectra behind such regions are therefore impor- 
tant. 

Since the ICM is confined in the proto-cluster poten- 
tial, the gas will move with the dark matter as the proto- 
cluster expands and collapses. We therefore do not expect to 
see wind-driven shells to form on the boundary of a proto- 
cluster. 



2.3 Connection with Lyman-Break galaxies 

Large numbers of star forming galaxies have been discov- 
ered at z J; 3 by colour selection using the Lyman-break 
technique (Steidel et al. 1996). The star-formation rates in 
these galaxies, as inferred from their UV-luminosities, are 
in the range 10 to IOOOM0 yr _1 (for a flat universe with 
Q, m — 0.3 and for a mean dust correction of ~ 7). The 
strong clustering of these galaxies is consistent with many 
of them being the progenitors of the early-type galaxies as 
seen in clusters today (e.g. Mo & Fukugita 1996; Governato 
et al 1998). Thus, if pre-heating of the proto-cluster medium 
occured at z <; 3, we may hope to see the kind of signatures 
discussed in the last subsection in the spectra of QSOs which 
lie behind a large cluster of LBGs. 

There is an interesting, although indirect, argument 
(Pagel 1999; see also Pettini 1999) that the metals returned 
from LBGs are mostly in hot gas. Indeed, if all the met- 
als ejected from LBGs were in cold gas, they should be 
observed as metal line systems associated with Lya ab- 
sorption systems. However, the observed total amount of 
metals in Lya absorption systems is much lower than that 
expected from the star formation inferred to occur in LBGs. 
A possible solution is that most metals are hidden in hot 
gas which would not produce significant Lya absorption. If 
this is the case, then metals are most likely contained in the 
proto-cluster medium, rather than being confined around 
individual galaxies, because radiative cooling would be very 
efficient if the gas were confined to such high density re- 
gions [see equation (H)] . Alternatively, the apparent scarcity 
of metals could be a selection effect, when metal-rich LBGs 
block background quasars due to dust obscuration. Clearly, 
observations of the proto-cluster medium around LBGs can 
help to distinguish between these two possibilities. 

If metal-enriched gas is being ejected from LBGs into 
the proto-cluster medium by winds, then we may hope to 
observe the signature of such gas outflows. Because of their 
high star-formation rate, LBGs are considered to be the 
high-redshift counterpart of local star burst galaxies (e.g. 
Heckman 2000). Therefore, we can use observed properties 
of local star bursts to predict the observable signatures of 
LBGs in the QSO absorption line spectra. As summarized in 
Heckman et al. (1999), local star bursts are regions of intense 
star-formation (more than IOM0 stars are being formed per 



year, in a region of ~ 1 kpc), which generally show gas out- 
flows with velocities in the range 400-800km s" 1 . These high 
outflow velocities do not strongly depend on the depth of the 
potential well (Heckman 2000) , and are sufficiently high that 
a galaxy can in principle pollute a large volume of the IGM 
in a Hubble time. The mass outflow rate is roughly propor- 
tional to the star- formation rate in the galaxy, M w oc M*, 
with the proportionality factor being typically a few. The 
wind consists mainly of cold clouds entrained in the out- 
flowing hot gas. The ratio between the mass outflow rate 
and the rate at which mass is being returned by supernova 
explosions and stellar winds is typically 5 to 10. However, 
only about 10 percent of the energy output from supernova 
explosions and stellar winds is in the kinetic energy of the 
cold flow; the rest is in the thermal and kinetic energy of 
the hot wind. 

We can construct a simple model (see also Wang 
1995ab) for the properties of the winds around LBGs, as- 
suming them to have similar outflow properties as star 
bursters. The density structure in a steady wind is 
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where u w is the wind velocity and u is the solid angle into 
which the wind blows. Inserting numerical values typical for 
a star burst wind gives for the hydrogen number density: 
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In the presence of an ionizing background, the neutral frac- 
tion of gas with clumping factor £ = (n^ n ) /(nHii) 2 is 
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where T4 ~ 1 is the temperature of the photo-ionized gas in 
units of 10 4 K, and the spectrum of the ionizing radiation is 
oc v~ x (e.g. Peebles 1993, section 23). Combining these two 
equations, we find for the neutral hydrogen column density 
at an impact parameter I 
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Thus, a QSO sight-line passing near a LBG will show a rel- 
atively strong, possibly broad Lya line, associated with the 
wind driven by the star burst. Complex velocity structure 
is to be expected if the outflow is clumpy. For an impact 
parameter smaller than ~ lOkpc, the QSO sight-line goes 
through optically thick clouds, and a Lyman limit system 
or a damped Lya system will be produced. 

The observational results by Pettini et al. (2000) for the 
spectrum of the gravitationally lensed galaxy MS1512-cB58 
(at z = 2.73) give direct evidence for such gas outflows with 
velocities ~ 200km s~ x in LBGs. However, common LBGs 
are too faint to yield spectra with high signal-to-noise ratios. 
Such outflows may be studied more efficiently is absorption 
towards background QSOs. 

Because of their large velocities, winds from LBGs may 
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drive shocks into the proto-cluster medium as they escape 
the potential wells of their parent halos. In the adiabatic 
phase, the shock front has a density contrast of about 4 
with respect to the ambient medium and a temperature T ~ 
fivl/2k ~ 2 x 10 e (u w /200km s _1 ) 2 K, where u is the mean 
molecular weight of the gas. Using equation (pi) we see that 
the shocked gas will remain very hot, as it is not be able to 
cool within a Hubble time. Thus, the winds may drive hot 
shells into the medium, which will eventually thermalize in 
the ambient proto-cluster gas. Since the gas in these shells is 
collisionally ionized, they are not expected to produce strong 
lines. 

In summary: a typical z ~ 3 quasar spectrum will in- 
tersect one proto-cluster per unit redshift. If the gas in this 
cluster is hot, as is expected from models that invoke pre- 
heating to explain the observed X-ray properties of clusters, 
then a large region with little Lya absorption would re- 
sult, which might also produce associated metal absorption 
in the form of Hubble broadened metal lines if the ICM is 
metal enriched. Such clearings in the forest would be appar- 
ent at sufficiently high redshifts, when the mean opacity in 
the Lya-forest is high. Proximity effects, due to AGN and/or 
galaxies, might make the clearing even more distinct. How- 
ever, if the proto-cluster gas is in a two-phase medium, with 
cold clouds pressure-confined in the hot medium, then these 
clouds may produce strong lines themselves. These could be 
distinguished from the more usual strong lines by the ab- 
sence of substructure, probed by higher-order Lyman lines, 
or by their small sizes, as measured from multiple quasar 
sight-lines. If the cold phase is metal enriched there would 
be a strong metal line associated with these lines. If proto- 
clusters usually contain a large number of Lyman-break 
galaxies, then QSO spectra taken through such a concen- 
tration can be used to test these predictions. Winds from 
LBGs are also expected to produce strong lines with col- 
umn densities > 10 15 cm -2 out to large impact parameters. 
These lines may have associated metal absorption if the wind 
is enriched. 



3 FEEDBACK FROM SMALL SYSTEMS 

Feedback on the scale of galaxy halos appears to be required 
for a variety of reasons. Without feedback, a major frac- 
tion of all baryons will cool and form stars in small halos 
at high redshift, exhausting the baryon supply to make up 
the IGM and the discs at lower redshift. This is known as 
the over cooling problem. The presence of an ionizing back- 
ground, for example due to AGN, may suppress the collapse 
of baryons into dark halos of circular speeds below ~ 20 — 30 
km s -1 (Efstathiou 1992, Quinn, Katz & Efstathiou 1996; 
Thoul & Weinberg 1996). For more massive systems, feed- 
back associated with star-formation is usually invoked to 
explain the low baryon content of dwarf galaxies and also to 
prevent forming galactic discs from loosing most of their an- 
gular momentum to their dark halo (Weil, Eke & Efstathiou 
1998). 

The basic picture for star-formation feedback is that hot 
bubbles blown by SNe and stellar winds overlap and power a 
galactic wind. Cosmic rays produced in the SNe shocks may 
also contribute (Breitschwerdt, McKenzie & Volk 1991). 
This feedback process heats cold gas which is then no longer 



available for star-formation, and may blow some fraction of 
the baryonic material out of the virial radius of the halo. In 
the next step of the merging hierarchy, the blown-out gas 
may once more be accreted. 

Semi-analytical models of galaxy formation (e.g. Kauf- 
man, White & Guiderdoni 1993, Lacey et al. 1993) require 
that a significant fraction of the total energy associated 
with SNe explosions be converted into thermal energy of the 
heated gas. However, numerical simulations of this process 
(e.g. Navarro & White 1994) fail to reach the required effi- 
ciency, because SN explosions invariably occur in dense gas 
which rapidly cools radiatively. Presently it is not yet clear 
whether the high efficiencies assumed in the semi-analytic 
models are unrealistic, or whether the simulations grossly 
underestimate the efficiency because they lack the resolu- 
tion to produce a multi-phase medium. 

Efstathiou (2000) recently proposed a model where 
feedback produces a multi-phase interstellar medium in 
which a low rate of star-formation can produce steady con- 
version of cold gas into hot gas, which then drives a wind. 
In this model, star formation is quenched because SNe ex- 
plosions limit the amount of cold gas available for star- 
formation. The model shows that a galaxy with circular 
speed ~ 50 km s _1 can expel 60-80 per cent of its gas over a 
time-scale of ~ 1 Gyr. Nath & Trentham (1997) also inves- 
tigated a model where an early population of dwarfs under- 
went strong mass loss through winds, and showed that this 
could pollute the low density IGM to a value close to those 
inferred from observations. The details of how such a wind 
interacts with the infalling gas may be very important for 
a realistic description of the feedback process, and also for 
estimating to what extent a forming galaxy can influence its 
surroundings. For example, Efstathiou (2000) argued that 
the wind might break-up into clouds which do not interact 
strongly with the accreting gas. Therefore, such a galaxy 
would undergo inflow and outflow simultaneously. 

In the next sections, we use our simulations to inves- 
tigate two very simple models for possible effects of galaxy 
feedback on the Lya-forest. In the first model we remove 
all the gas in the surroundings of galaxies so that it no 
longer contributes to the absorption. In this 'maximal feed- 
back' model, we assume that galactic outflows are very pow- 
erful and heat-up a significant fraction of the IGM surround- 
ing them. We will show that this has a dramatic effect on the 
number of strong absorption lines. In the next section, we 
examine the effects of galactic winds around dwarfs on the 
IGM. We show that shells formed from such a wind will over 
produce the number of Lya-lines, unless the wind breaks-up 
into clouds. 

3.1 Feedback and strong lines 

In the traditional picture of the Lya-forest, strong lines (col- 
umn density > 10 15 cm -2 ) are produced close to collapsing 
halos, which typically consist of a central concentration of 
high density, cooled gas which is surrounded by a hot halo 
of shocked gas. The hot halo itself is collisionally ionized 
and consequently contributes very little to the absorption. 
This is illustrated in Fig. |^, where we compare a stretch 
of absorption spectrum with and without the contribution 
from hot gas. Here, gas at a given density p is considered 
'hot' if its temperature is more than 20 per cent higher than 
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Figure 2. Simulated Lyct spectrum with (full line) and without (dashed line) the contribution from hot gas (see text for definition of 
'hot'). This stretch was carefully picked to show any difference at all between these two cases: hot gas contributes very little to the 
observed spectrum, as expected. 
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Figure 3. Simulated Lya spectra at z ss 3. Panel (a): simulated spectrum without feedback; panel (b): same as panel (a) but with 
feedback and scaled to have the same mean absorption as the spectrum in (a) (see text for details). Spectrum (b) contains many fewer 
strong lines. Panel (c) is the same as panel (b), but assuming the same amplitude of the ionizing background as in panel (a). 



the minimum temperature at that density. Even though this 
is a rather low limit, the spectrum shown in Fig. ^ is the 
only stretch of spectrum within a longer spectrum of length 
4 x 10 4 km s _1 (the typical length of a z = 3 QSO Lya spec- 
trum) where there is any noticeable contribution of the hot 
component. Consequently, as long as feedback on the hot 
component has no dynamical effects, it is unlikely to affect 
the observed properties of the Lya-forest. 

However, if feedback also affects the cold and dense 
gas in the surroundings of a halo, it may have an impor- 
tant effect on the Lya spectrum. We have run a standard 
friends-of-friends group finder (linking length 0.2 times the 
mean inter particle separation) on all stars in the simula- 
tion and identified 'galaxies' as groups of more than ten 
'stars' (one particle has mass 1.3 x 1O 6 M so the minimum 
stellar mass in a galaxy is 1.3 x 10 7 Mq). We then assumed 
that galaxy feedback is able to heat all the gas within a phys- 
ical radius R around the galaxy to such high temperatures 
that it no longer contributes to the absorption. We choose 
R = 150kpc(i.e., a significant fraction of the virial radius 



of a large galaxy), in which case the volume filling factor, 
n ga i(47r/3)7i 3 = 1, where n ga i is the number density of galax- 
ies. (The true volume filling factor of hot gas is then much 
lower, since these clustered halos overlap considerably.) 

Fig. |^ compares a stretch of simulated spectrum (at 
z ~ 3, panel a), with a simulated spectrum with feedback 
(panel b). The simulated spectra are scaled to have the same 
mean absorption as the spectrum of QSO 1422+231 (Rauch 
et al. 1997). Feedback significantly decreases the number of 
strong lines, since the gas that originally produced these 
lines is now assumed to be hot (compare e.g. the lines at 
v « 2000 km s" 1 in panels a and b). Conversely, weak lines 
in panel (a) appear much stronger in panel (b), since the 
spectra are scaled to the same mean absorption. Spectrum 
(a), without feedback, looks similar to the observed spec- 
trum of QSO 1422 in terms of the number of strong lines, 
whereas spectrum (b), which mimics strong feedback, has 
far fewer strong lines. Panel (c) is a spectrum with feedback, 
but not scaled to the same mean absorption, i.e., assuming 
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Figure 4. Comparison of the cumulative fraction of pixels with 
flux greater than some value. Full line: spectrum of QSO 1422, 
with la error contours (dashed lines); filled circles: simulated 
spectrum without feedback; squares: simulated spectra with feed- 
back due to dwarf galaxies (filled squares) and due to all galaxies 
(open squares), see text for details. Feedback has a significant 
effect on the number of pixels with high absorption, flux < 0.2. 



the same amplitude for the ionizing background as in panel 
(a). This illustrates the effect of the rescaling. 

Indeed, the removal of saturated regions due to feedback 
has a clear observational signature in a more detailed anal- 
ysis. Fig. ^ compares the one-point functions (histogram of 
the flux per pixel) of the spectra shown earlier. Whereas 
the simulated spectrum without feedback reproduces the 
observed one-point function very well, the simulations that 
mimic feedback produce significantly fewer pixels with high 
absorption (flux < 0.2). Two such simulations are illustrated 
in Fig^, the one labelled 'dwarfs' includes feedback from 
galaxies with mass 1.3 x 10 7 < M/M Q < 1.3 x 10 8 , the 
one labelled 'all' includes feedback from all galaxies above 
1.3xl0 7 M o . In both cases, the radius of the sphere R around 
the galaxy from which the cold gas is removed, prior to com- 
puting the mock absorption spectra, equals R = 150kpc, and 
the volume filling factors are then w 0.1 and ~ 1, for dwarfs 
and all galaxies, respectively. [We recall that these volume 
filling factors are simply the ratio of the sum of the vol- 
umes of all bubbles to the total volume. The true fractions 
of volume encompassed by the bubbles, i.e. taking into ac- 
count that the bubbles overlap, are 6.5 and 30 per cent, for 
the dwarfs and all galaxies respectively. If we increase R to 
220kpc for the dwarfs simulation - in which case the vol- 
ume filling factor increases to ~ 1 (true volume filling factor 
38 per cent), the one-point function (as well as the number 
of strong lines, to be discussed next), hardly changes.] Cold 
gas surrounding the more massive galaxies produces a rela- 
tively large fraction of the saturated pixels, which explains 




lo gio N m (cm- 8 ) 

Figure 5. Cumulative number of Voigt-profile lines per unit ve- 
locity (km s ), with column density larger than some value. 
Symbols are the same as in Fig. ^. The simulations that mimic 
feedback (open and filled squares) produce significantly fewer 
lines than the simulation without feedback (filled circles). The 
latter reproduces the observations (full line) very well. Dashed 
lines indicate the 25 and 75 per centiles for the observations, de- 
termined using bootstrap resampling. 

the difference between imposing feedback on all galaxies, or 
only on the small ones. 

Cen (1997) argued that the high-end of the Lya optical 
depth distribution is quite sensitive to the amplitude of the 
power spectrum, because it probes the tail of the (e.g. Gaus- 
sian) density probability distribution. If feedback effects are 
strong then the ability of this statistic to distinguish between 
cosmological models may be partly compromised. 

The consequence of having fewer strong lines, in the 
presence of feedback, can also be seen when plotting the 
column density distribution, d 2 N/dNni/dV, the fraction of 
lines per unit column density, per unit spectral length in 
km s~ . We have fitted Voigt profiles to both simulations 
and to the observed spectrum of QSO 1422, using the same 
fitting algorithm (an automated version of VPFIT; Webb 
1987, Carswell et al. 1987). Fig. ^ shows that the simula- 
tion without feedback reproduces the observed numbers of 
strong lines very well, whereas the simulations that mimic 
feedback significantly under produce the number of strong 
lines. Note that the shape of the column-density distribution 
is sensitive to the amount of large-scale power (e.g. Gnedin 
1998). This could be used to distinguish the effects of feed- 
back from those of large-scale power, on the column-density 
distribution. 

The above considerations show that this type of max- 
imal feedback, occurring in a spherical region around the 
forming halo can be ruled out. However, the distribution 
of gas that gives rise to absorption lines is actually far from 
spherical. In fig. 0we plot the matter distribution of a system 
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Figure 6. Geometry of the structure that produces a strong Lya line. Top left panel: Lyo absorption spectrum (flux versus velocity) 
of a line with column density 1.4 X 10 16 cm -2 (full line) and the corresponding Ly/3 profile (dashed). Panels a, b and d show the gas 
distribution in a thin slice along the sight-line (the dashed horizontal line in panels a-e) that produces this strong line. Plotted are gas 
particles with p > (p) and temperature T < 2.5 X 10 4 K (panels a-b) and hot gas particles with T > 2.5 X 10 4 K (panel d). Positions 
are labeled in Hubble velocity v = H(z)r. The slice in the y-direction (perpendicular to the page) is 10 km s _1 (20 km s _1 for panel 
d). Star particles are plotted as stars in panel (c). In panel (a), v z refers to the rcdshift space velocity, Hr + v pec - Panel (e) plots the 
corresponding dark matter particles. The absorption line is a superposition of three to four filaments which overlap in redshift space 
due to the high peculiar velocities induced by the nearby large mass concentration. These high velocities make it difficult to identify the 
structure in redshift space (panel a) with the corresponding filaments in real space (panel b). The massive system at the top is embedded 
in hot, shocked gas, but there is very little hot gas associated with the filaments that produce the line. Note that the gas density is an 
interpolation over the plotted SPH particles, so that particles not directly hit by the plotted line-of-sight may still contribute to the 
absorption. 



that produces a strong line with column density 1.4 x 10 
cm~ 2 . This line is a superposition of the absorption from 
three to four different filaments, which blend together in 
redshift space due to peculiar velocity gradients. The dis- 
tribution of the cold gas which is responsible for the line is 
filamentary. Presumably, a galactic outflow from the stars 
near the dense knot at v z ~ 600km s _1 , which contributes 
most to the absorption, would preferentially occur perpen- 
dicular to the filament along the path of least resistance, in 
which case the line may not suffer strongly from feedback 
after all. But the consequence of directing all the feedback 
momentum into the low-density surroundings, is that the 
accretion rate onto the halo, which is dominated by infall 
along the filament, cannot be strongly affected by the feed- 
back. That is, if feedback is able to limit the accretion rate 
onto dwarf halos, then it will necessarily also influence the 
properties of strong Lya absorption lines. Note that the in- 



dividual components that give rise to the strong line can be 
recognized in the higher-order Ly/3 line, if it is not saturated. 
Therefore, higher order Lyman-lines should be good probes 
of the state of the cold gas in the surroundings of forming 
halos, and can be used to examine whether feedback is able 
to significantly stir the accreting gas. 

Clearly, a more detailed analysis, for example based on 
probing substructure in strong lines using higher order tran- 
sitions, is required, but it does seem safe to conclude that the 
properties of strong lines may be a valuable probe of feed- 
back round forming galaxies. A detailed comparison of sim- 
ulations without feedback with high-resolution observations 
may shed light on whether galaxy feedback is an important 
factor in determining the properties of strong lines. 
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Figure 7. Dynamics of a wind with mass-loss rate M w = 
O.5M0/yr and wind-speed 100km s — 1 , emanating from a galaxy 
with circular velocity 50km s _1 into the IGM with sound speed 
20km s — x at a redshift z = 3. (We assumed in addition h = 0.65, 
Qt = 0.05 and Q m = 0.3 for the cosmological parameters.) Panel 
(a) depicts the velocity of the shell (in units of the wind velocity) 
versus its position (in units of the virial radius R200 = 17.5kpc of 
the halo). Panel (b) shows the mass of the shell, in units of the 
initial mass. The dashed curves neglect the forces acting on the 
shell. Panel (c) shows the magnitude of the dimensionless forces 
acting on the shell. From top to bottom, the curves denote the 
total force (full line), gravity (short dashed), ram-pressure force 
of the surrounding IGM (long dashed), ram-pressure of the wind 
(dot-short dashed), and the thermal pressure of the IGM (short 
dashed-long dashed. 



3.2 Winds from dwarf galaxies 

In the previous section we discussed the effects of feedback 
on the immediate surroundings of the galaxy. We showed 
that feedback cannot influence all the cold gas around form- 
ing halos out to ~ 150kpc without significantly affecting the 
properties of strong lines as well. Here we discuss in more 
detail whether outflows from dwarf galaxies may themselves 
be responsible for producing absorption lines. 

The mass-to-light ratio in dwarf galaxies is typically of 
the order M/L ~ IOOMq/Lq (e.g., Carignan & Freeman 
1988, Lake, Schommer & van Gorkom 1990), much higher 
than in spirals or ellipticals. Presumably, a wind associated 
with feedback from star-formation was able to blow most 
of the gas out of the shallow potential well of the galaxy. 
Indeed, in the model of Efstathiou (2000) , a galactic wind is 
able to remove 60-80 per cent of the gas contained in a dwarf 
galaxy, over a time-scale of order lGyr. In these models, the 
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Figure 8. Stagnation radius of the shell (in units of the virial 
radius i?200 of the halo) versus circular velocity of the halo, for 
wind velocities of 50, 100, 150 and 200km s -1 (panel a), and 1, 2, 
3, 4 



outflow velocity at the virial radius is typically of order a 
few times the escape velocity, which in turn is typically a 
few times the circular speed (« C irc) of the halo. In this case, 
the wind may evacuate a cavity in the IGM, surrounded by 
a shell. 

The momentum conservation equation for a spherical 
shell with mass M s and velocity V a , which is being blown 
into the surrounding IGM by a galactic wind with mass-loss 
rate M w and velocity V m is 

j t (M s V s )= M W V W (1- ^ L )-4irR 2 3 p{C 2 /~ f + V{V S + V)} 



GMg M s 

m 



(21) 



where R s is the current position of the shell. The first term 
is the momentum accreted from the wind. The second term 
is the deceleration of the shell due to the pressure of the sur- 
rounding IGM (density p(R a ), sound speed C and adiabatic 
index 7), and due to the sweeping-up of matter which rains 
onto the expanding shell with velocity V. Finally, the last 
term is the gravitational deceleration due to the mass M g 
interior to the shell (G denotes Newton's constant of grav- 
ity). The mass accretion rate of the shell due to accretion 
from the wind and the sweeping-up of the IGM is 



d_ 



M a = Af«(l - + 4tyR 2 s P (Vs + V) . 



(22) 



The initial conditions can be set at the virial radius -R200 
of the halo, defined as usual as the radius within which the 
mean overdensity is 200 times the critical density. Mo, Mao 
& White (1998) give the relation between the virial radius, 
the enclosed mass M200, the circular velocity V c i rc and the 
Hubble constant at redshift z as 

Rwo = wWli) ' (23) 
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We will assume V 3 = V w at R a — R200 and denote the 
fraction of (baryonic) mass swept-up by the wind by the 
time it reaches the virial radius by a < 1. We also assume 
the wind-speed to be constant. Finally we need to specify 
the density profile of the material surrounding the galaxy, 
for which we assume 



p{R) = 



200 / -R200 , 
3 V R 1 



(25) 



where p c = 3H 2 /8tvG is the critical density. The gas density 
is Q b /Q m p. These equations are valid at high redshift where 

(p)~Po- 

Equations ( J2 1| ~ |22| ) can be cast in dimensionless form, 
using -R200, V w and M200 as the units of distance, velocity 
and mass, respectively. Using lower-case letters to denote 
dimensionless variables we find 
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(m s v s ) = 



(1 - Vs) 



M200 

tt b r 2 s {r- 2 + 3/200) [(v + v 3 )v + c 2 / 7 ] 



; (1 + x ~ m s )m s /r s 



(26) 



dr 



MwU 

M200 



(1 - O + n b r 2 s (r7 2 + 3/200) (v + v s ) 

(27) 



where x = ( r s — 1) + (ff — l)/200 and the time variable 
t = t/t* with t* = R2oo/V w . 

Fig. illustrates the dynamics of a wind emanating 
from a small galaxy with circular velocity 50km s . We 
assumed that feedback is very efficient and the initial shell 
mass M s = a(fi;,/f2 m )M2oo with a = 1. We further as- 
sumed that gas is accreting onto the galaxy with velocity 
V — Vcirc and integrated the equations until the stagna- 
tion point where the velocity becomes sonic, V s = C. The 
deceleration of the shell is dominated by gravity (mostly 
from the halo) in this case, with a non-negligible contri- 
bution from the ram-pressure of the infalling gas. Nath & 
Trentham (1997) followed Tegmark, Silk & Evrard (1993) 
and neglected the gravity of the halo, as did Ferrara, Pettini 
& Shchekinov (2000) when they estimated the distance to 
which a wind-blown shell will travel and pollute the IGM. 
In our case, the shell travels out to ~ 1.8 times the viral ra- 
dius, thereby doubling its mass. The dashed lines in panels 
(a) and (b) illustrate the case where forces are neglected, in 
which case the shell reaches a considerably larger distance, 
sweeping-up ~ V w /C times its initial mass. 

Fig. ^ illustrates the stagnation point of the shell for 
different wind velocities, as a function of the circular velocity 
of the halo. We have assumed that a wind with mass-loss 
rate M m = O.5M0 yr _1 blows for 0.5Gyr (so that it can in 
principle remove a substantial fractrion of the baryons over a 
long time; see also Efstathiou 2000) , and plot the position of 
the shell when it stagnates, V s = C, or until t = l.OGyr. In 
small halos, winds of 100 — 150 km s _1 are able to expell gas 
up to 3-6 times the virial radius, but the increased gravity in 
more massive halos prevents the shell from traveling far. On 
the other hand, if the wind velocity is a given fraction of the 
circular speed, then slow winds blow further in more massive 



halos, since such winds can blow for longer before the shell 
speed becomes equal to the sound speed of the IGM. 

The virial radius -R200 of a halo is related to its circular 
velocity by (Eqs. 00) 



-R200 = 17kpc ■ 



Kir 



_f^Y 1/2 fl+±y 3/2 
1 V0.3J V 4 J 



50km s 

and the corresponding virial mass 

v- \ 3 /n \~ 1 ' 2 
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(29) 



Therefore for the wind parameters in Fig. (|7j) the final shell 
radius R s ~ 30kpc and the mass in the shell M s ~ 3.3 x 
1O 9 M , where we assumed that Q, b /Q. m = 0.05/0.3 = 0.17. 

The cross-section of shells of such a size is quite high, 
with a mean number of intersections per unit redshift of 



dN 
dz 



50 



( Rs 



0.5(/i" 1 Mpc) 3 V 30kpc 



f— r 

V0.65/ 



V03 



-1/2 



1/2 



(30) 



where no ~ 0.5(/i -1 Mpc) -3 is the co-moving number- 
density of dwarfs (as determined from our galaxy-cluster 
normalized simulation). If such a shell produces absorption 
lines, then a significant fraction of all Lya lines could be due 
to such shells. (Az = 3 quasar spectrum contains typically 
of order 300 lines above a column density of 10 13 cm -2 per 
unit redshift). 

Note that one would expect such wind blown bubbles to 
produce two lines, when punctured by a sight-line, separated 
by ~ 100km s _1 . Cristiani et al (1997) analyzed Ly a- lines 
of several high-resolution quasar spectra and concluded that 
the two-point correlation function of the lines shows clus- 
tering at such scales. In contrast, Pettini et al (1990) and 
Stengler-Larrea & Webb (1993) claimed a non-detection of 
such a signal. Part of the confusion might be caused by the 
non- uniqueness of Voigt-profile fits, especially in the pres- 
ence of severe line-blending occuring at high redshifts. A 
uniform comparison of high-resolution spectra with simula- 
tions is required to check for the reality of such a signal. 

Below we will estimate the column density of such a 
shell, but first we would like to point out that there are 
significant uncertainties in the estimate of the shell's radius. 
Dwarf galaxies in simulations are invariably embedded in fil- 
aments at over densities of a few, which have significant non- 
spherical infall velocities which will oppose outflows. There- 
fore a more realistic description of the shell's motion should 
take into account such a filamentary matter distribution. 

More importantly, it is not clear how efficiently the out- 
flowing gas will couple to the surrounding IGM. The wind 
may break-up into clouds which coast to large distances 
without sweeping-up a significant fraction of the ambient 
IGM. In that case, although some of the outflowing gas may 
still reach large distances, the swept-up mass may be much 
less than that estimated above. If we assume that only a frac- 
tion / of the surrounding IGM is swept-up, with the rest in 
cold clouds, then the total column density of the remaining 
gas 
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N H ~ 3.5 x 10 la / 



M s ( R s V 
3 x 10 9 M o I 30kpc J 



(31) 



Since the wind velocity (~ 100km s _1 ) is much larger than 
the sound speed of the IGM, the wind results in a strong 
shock as it runs into its surroundings. In the adiabatic phase, 
the density of the shocked shell is then about 4 times that of 
the ambient medium (ni), while the temperature is about 
T 2 ~ 5 x 10 5 (u w /200km s _1 ) 2 K. Gas of such temperature 
and density, and with a metallicity Z > O.OIZr, has a cool- 
ing time shorter than its age [see equation (pi). At these 
temperatures, cooling will be dominated by metals, hence 
A23 3> 1.]. Thus, as the shock propagates into the medium 
and sweeps-up gas, it can cool down to ~ 10 4 K and become 
dense, with a density 713 ~ Mfni ~ lOOni (~ 10~ 3 cm -3 at 
z ~ 3), where A4i is the initial Mach number of the shock. 
Under these assumptions, the column density of the shell is 



Nm 



1.4 x 10 16 / 



J21 



( n * ) 
Vl0^ 3 cm- 3 / 



" 3.5 x 10 19 cm- 2 ' ^ 

where we used equation ([[9]) for converting the total to the 
neutral column density. 

The physical parameters for the winds from a dwarf 
versus those from a LBG, discussed in Section 2.3, are quite 
different. For the assumed mass-loss rates (0.5 for the dwarf 
versus lOOM© yr _1 for the LBG) and wind velocities (100 
versus 200km s _1 ) for these systems, the wind densities dif- 
fer by a factor 1/100. Consequently for the LBG, it is the 
dense, cold wind that produces most of the absorption, since 
the swept-up hot shell cannot cool. The reverse is true for 
the wind from the dwarf galaxy, for which the tenuous wind 
produces no significant absorption, but the dense swept-up 
shell produces lines. 

If the wind interacts strongly with the IGM, / ~ 1, 
the high column density of the lines given by equation ( |32| ) 
combined with the rate of incidence of the shells, dN/dz ~ 
50, predicts of order 100 lines with columns larger than 
10 16 cm -2 per quasar spectrum, which is already much 
more than the typical total number ~ 20 of such lines in 
an observed spectrum. Therefore, these considerations re- 
quire / < 1, that is the wind cannot retain its integrity 
but must break-up into clouds. These clouds, which may 
have higher column-densities, are then not subject to the 
momentum conservation Eq. ( [ill ) an( l can coast to large 
distances r ~ 0.3 Mpc if they move at 100km s _1 over 
the Hubble time, Eq. ([s]). Since these clouds are polluted 
by metals of the SNe explosions that drive them out of the 
galaxy in the first place, such winds should be able to de- 
posit heavy elements throughout the IGM. Clearly a much 
more detailed model is required to compute the covering 
factor of such clouds. Alternatively, the matter density in 
the surroundings of the halo may limit the expansion of the 
shell much more than we have computed here, for exam- 
ple because the infall velocities are higher than we have as- 
sumed, or because the distribution is non-spherical and the 
winds are tunneled into the low-density surroundings with- 
out sweeping-up much matter. 

Finally we would like to make a remark on the effects of 
such winds on neighbouring halos. In numerical simulations 
of hierarchical clustering, each dark matter halo always con- 



tains a swarm of smaller halos down to the resolution scale 
(e.g. Moore et al 1999). A galactic wind would presumably 
remove most of the gas from such sub-clumps, leaving these 
sub-structures mostly dark. 

In summary: strong Lya lines may be excellent probes 
of feedback effects on galactic scales. If galactic feedback is 
assumed to heat a large fraction of the cold gas in the sur- 
roundings of forming halos, then many fewer strong lines are 
produced than are observed in QSO spectra. Shells formed 
from gas swept-up by galactic winds from dwarf galaxies 
would over-produce the number of strong absorption lines. 
Therefore we argued that these winds should break-up into 
clouds, which can then coast to large distances and pollute 
the IGM with metals, or be funneled into the low-density 
surroundings without producing a strong line. 



4 SUMMARY 

The paradigm in which the Lya-forest is produced by 
the low-density photo-ionized IGM, which traces the struc- 
ture of the underlying dark-matter perturbations, has been 
very successful in explaining properties of observed high- 
resolution quasar spectra. However, models of galaxy and 
galaxy cluster formation require large amounts of feedback 
in order to match the observations. We have investigated 
possible observational signatures of this type of feedback in 
Lya spectra. 

Pre-heating is often invoked to explain the observed X- 
ray properties of galaxy clusters. A QSO sight-line inter- 
sects of order one such proto-cluster per unit redshift. If 
the intra-cluster gas is hot, it will produce a large gap with 
little absorption in the spectrum, which could be detected 
at high redshift, z ~ 3 — 4, when the Lya-forest opacity 
is high. If the proto-cluster gas is in a two-phase medium, 
the gas in pressure-confined cold clouds could produce a sin- 
gle strong line with associated metals. Such a model can be 
tested observationally, by probing the sub-structure of such 
absorbers using e.g. higher-order Lyman lines, or by using 
multiple sight-lines. Dense groups of Lyman-break galaxies 
at z ~ 3 — 4 are likely sites of proto-clusters, and QSO 
spectra taken through such a structure could be used to 
put limits on cluster pre-heating. Winds associated with 
star formation in Lyman-break galaxies will also produce 
lines with column-densities > 10 15 cm -2 , plausibly associ- 
ated with metals escaping from the galaxy. Such lines can 
be used to study the feedback process in these galaxies. 

Feedback on galactic scale will influence the properties 
of strong Nm ~ 10 16 cm -2 Lya lines, since these are pro- 
duced close to forming halos. If feedback is able to heat a 
significant fraction of the cold gas in the surroundings of a 
galaxy, then many fewer strong lines will be produced. How- 
ever, the outflows may be preferentially directed into the 
low-density surroundings, in which case strong lines may 
survive. Shells associated with winds from dwarf galaxies 
will over produce the number of absorption lines, unless the 
wind breaks-up into clouds. These clouds can coast to large 
distances and pollute the IGM with metals. Alternatively, 
the matter surrounding the galaxy may funnel the wind into 
the low-density IGM, without producing lstrong lines. 

These simple arguments suggest that feedback may 
have an observable effect on the Lya-forest. Some of our 
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predictions on cluster feedback can be tested by correlat- 
ing Lya-absorption spectra with the Lyman-break galaxies 
along the line-of-sight. Simulations of dwarf galaxy forma- 
tion that include a wind, can be used to examine the inter- 
action of a galactic wind with accreting gas. 
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